Collective modes and sound propagation in a p-wave superconductor: Sr2Ru04 
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There are five distinct collective modes in the recently discovered p-wave superconductor Sr2Ru04; 
phase and amplitude modes of the order parameter, clapping mode (real and imaginary), and spin 
wave. The first two modes also exist in the ordinary s-wave superconductors, while the clapping 
mode with the energy \/2A(r) is unique to Sr2Ru04 and couples to the sound wave. Here we 
report a theoretical study of the sound propagation in a two dimensional p-wave superconductor. We 
identified the clapping mode and study its effects on the longitudinal and transverse sound velocities 
in the superconducting state. In contrast to the case of '^He, there is no resonance absorption 
f-H I associated with the collective mode, since in metals tt;/(t;ir|q|) <^ 1, where vp is the Fermi velocity, 

Q ■ q is the wave vector, and u is the frequency of the sound wave. However, the velocity change in the 

O ' collisionless limit gets modified by the contribution from the coupling to the clapping mode. We 

compute this contribution and comment on the visibility of the effect. In the diffusive limit, the 
^ l' contribution from the collective mode turns out to be negligible. The behaviors of the sound velocity 

^ ' change and the attenuation coefficient near Tc in the diffusive limit are calculated and compared 

^ , with the existing experimental data wherever it is possible. We also present the results for the 

' attenuation coefficients in both of the collisionless and diffusive limits at finite temperatures. 
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PACS numbers: 74.20.-z, 74.25.Ld, 74.25.-q 



I. INTRODUCTION 



' Shortly after the discovery of superconductivity in Sr2Ru04, the possibility of spin triplet pairing was discussed. 

^ [Ql Possible pairing symmetries were also classified based on the crystal symmetry. |^ On the experimental front, 

T-H ' there have been attempts to single out the right pairing symmetry among these possibilities. Recent measurement of 

CO , ^^0-Knight shift in NMR for the magnetic field parallel to the a — b plane showed no change across which can be 

] • taken as the evidence of the spin triplet pairing with c?-vector parallel to the c-axis. |^ Here d is called the spin vector 

^ , which is perpendicular to the direction of the spin associated with the condensed pair Q. /iSR experiment found 
spontaneous magnetic field in the superconducting Sr2Ru04, which seems to indicate broken time reversal symmetry 

^\ , in the superconducting state. ||^ These experiment may be compatible with the following order parameter j^] 

■ A{k) ^ Ad{ki±ik2), (1) 



I ■ where A is the magnitude of the superconducting order parameter. Notice that this state is analogous to the A phase 
' O \ of ^He and there is a full gap on the Fermi surface. 

On the other hand, there also exist experiments that cannot be explained by a naive application of the order 
parameter given by Eq.(l). Earlier specific heat measurement found residual density of states at low temperatures 
below Tc ^ , which provokes the ideas of orbital dependent superconductivity |^] and even non- unitary superconducting 
state However, more recent specific heat experiment on a cleaner sample reports no residual density of states 
and it was found that the specific heat behaves as at low temperatures. P] This result stimulated a speculation 



?H ' about different order parameters with line node |10|| . However, since there are three bands labeled by a, (3, and 7 
. . - which cross the Fermi surface, it is not yet clear whether the order parameter given by Eq.(l) is compatible with more 
recent specific heat data or not. For example, it is possible that the pairing symmetry associated with the 7 band is 
still given by Eq.(l) while the order parameter symmetry associated with a and (3 bands can be quite different. In 
this case, the low temperature specific heat will be dominated by the excitations from a and /3 bands. In order to 
resolve the issue, it is important to examine other predictions of the given order parameter and compare the results 
with future experiments. 

One way of identifying the correct order parameter among possible candidates is to investigate the unique collective 
modes supported by the ground state with a given pairing symmetry. The observation of the effects of these collective 
modes would provide convincing evidence for a particular order parameter symmetry. If we assume that the order 
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parameter of Eq.(l) is realized in Sr2Ru04, the superconducting state would support unique collective modes, the 
so-called clapping mode and spin waves as well as the phase and amplitude modes of the order parameter which 
exist also in s-wave superconductors. Previously we studied the dynamics of spin waves QJl^. A possible way to 
distinguish the order parameter of the 7 band from those of the a (or f3) band was also proposed in the context of 
spin wave dynamics. 

In this paper, we study the dynamics of the sound wave and its coupling to the clapping modes assuming that the 
order parameter is given by Eq.(l). As in ^He, only the clapping mode can couple to the sound wave and affects its 
dynamics. Here we study the sound velocities and attenuation coefficients of the longitudinal and transverse sound 
waves. In particular, we identify the clapping mode with the frequency, uj — \/2A(T), and examine the effects of this 
mode and disorder on the sound wave propagation. 

In a recent paper, Higashitani and Nagai obtained the clapping mode with the frequency, u) = v^A, and 
discussed the possible coupling to the sound wave independent of us. It is, however, important to realize that the 
coupling to the sound wave is extermely small because C/vp 1 in metals, where C is the sound velocity. Indeed 
the recent measurement of the sound velocity in the normal and superconducting states of Sr2Ru04 reported in 
shows that C/vp <C 1. They measured the sound velocities of the longitudinal modes, Cn (q, u || [100]) and C33 
(q,u II [001]), and the transverse modes, C44 (q || [100], u || [001]) and Cee (q || [100], u || [010]), where q and u 
are the directions of propagation and the polarization of ultrasound, respectively. They found that the longitudinal 
sound velocities, Cn and C33, decrease with a kink at T = Tc, while the transverse sound velocities do not exhibit 
any effect of the onset of superconductivity. We estimate from their experimental data that Ci/vp ~ 10~^, where C; 
is the longitudinal sound velocity. It can be also seen that the transverse sound velocity, Ct , is much smaller than the 
longitudinal one | |T^ ]. 

Incorporating the correct limit, C/vp <C 1, we obtained the sound velocities and attenuation coefficients for both 
coUisionless and diffusive limits. In the diffusive limit, the quasi-particle scattering due to impurities should be 
properly taken into account. One can show that, in a metal like Sr2Ru04, the coUisionless limit is rather difficult 
to reach because it can be realized only for uj ~ 0(1) GHz. For more practical range of frequencies, kHz — MHz, 
the diffusive limit may be easier to achieve. On the other hand, we found that it is much easier to see the effects of 
the coupling between the sound waves and the clapping mode in the coUisionless limit. Therefore it is worthwhile to 
study both regimes. 

Here we summarize our main results. 

A. CoUisionless limit 

In the absence of the coupling to the clapping mode, the longitudinal sound velocity decreases in the supercon- 
ducting state because the effect of the screening of the Coulomb potential increases, which happens in the s-wave 
superconductors as well. However, one of the important features of the p-wave order parameter in consideration is 
that the sound wave can now couple to the clapping mode. This effect is absent in s-wave superconductors. One can 
show that, among longitudinal waves, Cn mode can couple to the clapping mode, but C33 mode cannot. We found 
that the longitudinal sound velocity Cn decreases as 



2 \vp J ' ■' A{l + {2A{T)/vpqy}' 



(2) 



where A; is the couping constant and / is the superfluid density. 6C1/C1 is the relative shift in the sound velocity. 
We estimated the frequency regime where one can observe the effect of the clapping mode and found that the effect 
is visible if vp\q\ ~ 2 — 3A(0). This implies that uj ^ 0{1) GHz. Since C33 does not couple to the clapping mode, 
the velocity change is simply given by 



SCf^ _ ^33 



(3) 



Since the velocity of the transverse wave is much smaller than that of the longitudinal one and the coupling to 
the electron system is weaker than the longitudinal case as well, we expect that the change of the transverse sound 
velocity is hard to observe. In order to complete the discussion, we also present the results for these small changes in 
transverse velocities. Here only Cee mode couples to the clapping mode, and C44 mode does not. We found 
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2 \vp J \ ■' 4{1 + {2A/vpq)^ 



(4) 
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(5) 



where the At is the transverse couphng constant. 

We also found that the leading contribution to the attenuation coefficient is the same as that of s-wave supercon- 
ductors in the coUisionless limit. 
B. Diffusive limit 

This case corresponds to Ui?|q| <C T, where T is the scattering rate due to impurities. As in the case of the 
coUisionless limit, in principle Cn and Cge modes couple to the clapping mode, but it turns out that the effect is almost 
impossible to detect. Neglecting the coupling to the clapping mode and working in the limit 'iirTc ^ 2r ^ i;i?|q|, w, 
we obtain the following results near Tc. 



47r3 T / T 

7C(3)r \~t'. 
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7C(3)r 
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where a and Q!„ are the attenuation coefhcients in the superconducting and the normal states respectively. The 
shift of the sound velocity and attenuation coefficient decrease linearly in (1 — T /Tc) as T ^ Tc. This result for the 
longitudinal sound wave is consistent with the experimental observation reported in [p^ 

In the case of the transverse sound waves, the leading behaviors of the sound velocity and the attenuation coeffient 
can be obtained simply by replacing Ci and A; by Ct and At in the diffusive limit. However, the absolute value of the 
transverse sound velocity is much smaller than the longitudinal one and the coupling to the electron system is also 
much weaker than the case of the longitudinal sound waves. Thus, it would be hard to observe any change at T — Tc 
for the transverse wave. This may explain the experimental finding that the transverse velocity does not show any 
change across Tc. We also obtained the attenuation coefRcient for all temperatures below Tc. It is given by Eq.(|3^) 
and Fig. 2 shows its behavior. 

The rest of the paper is organized as follows. In section II, the clapping mode is briefly discussed. In section III, 
we provide a brief summary of the formalism used in Ref. p^ ] to explain how the sound velocity and the attenuation 
coefficients are related to the autocorrelation functions of the stress tensor. We present the results of the study on the 
sound propagation in the coUisionless and diffusive limits in sections IV and V, respectively. We conclude in section 
VI. Further details which are not presented in the main text are relegated to the Appendix A and B. 



II. COLLECTIVE MODES IN SR2RUO4 

As in the s-wave superconductors, the phase and amplitude modes of the order parameter also exist in the p-wave 
superconductors. On the other hand, due to the internal structure of the Cooper pair in the p-wave superconductor, 
there exist other types of collective mode associated with the order parameter. The nature of these modes is determined 
by the structure of the order parameter. 

There are collective modes associated with the oscillation of the spin vector d, which we have already discussed 
in ||ll| and There exists another collective mode associated with the orbital part. Using the notation e^"^ = 

(fci ±ifc2)/|k|, the oscillation of the orbital part e**"^ e^*"^ gives rise to the clapping mode with to = -\/2A(T). This 
mode couples to the sound waves as we will see in the next section. Therefore, the detection of the clapping mode 
will provide a unique evidence for the p-wave superconducting order parameter. The derivation of the clapping mode 
and the coupling to the sound wave is discussed in Appendix A. 



III. DYNAMICS OF SOUND WAVE VIA STRESS TENSOR 



In ordinary liquids, the sound wave is a density wave. In superconductors, the density is not only coupled to the 
longitudinal component of the normal velocity, but also to the superfluid velocity and to temperature. The role of 
these couplings and their consequences in the dynamics of sound wave can be studied by looking at the autocorrelation 
function, ([t^ , Ty]), of stress tensor, t^ . 

([t,j, r,,])(r - r', t - t') = -ie{t - t'){[n,{v, t),n,{r', t')]) , (7) 
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where 



(v-v)av-Vy,- ^,t 



2i 



2im 



(8) 



Here tp'l is the electron creation operator with spin a. 

The other operators whose correlation functions are needed for the ultrasonic attenuation and the sound velocity 
change are the density operator 



(9) 



and the current operator 



2im 



V'J(r,t)V'a(r',i) 



(10) 



Assuming that the wave vector of the sound wave is in the x direction, q = qx, the sound velocity shift, SC, at low 
frequencies can be computed from 



SQ ^ Ci{uj)~Ci 
Ci Ci 

6a ^ Ct{u;)-Ct 
Ct 



Ct 



<^=Ci\q\ 



^=Ct|q| 
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■Re([/i;,/i;])(q,w) 



■Re{[ht,ht])iq,Lo) 



a)=Ct|q| 



(11) 



where 



hi[v,t) = -Ta;:r(q, i) n{v,t), 

(jj q 

ht{r,t) = -T^y{(i,t) - ■mjy{r,t), 



(12) 



Here C/ and Ct represent the longitudinal and transverse sound velocities in the normal state respectively, mjon and 
771 are the mass of ions and the mass of electron, respectively. On the other hand, the attenuation coefficient, a, at 
low frequencies is obtained from 



ai 



niionCl 



lm{[hi,hi]){q,uj) 



OJ 



at 



=C,|q| 



rriionCt 



Im([/it,/i(])(q,w) 



(13) 



=Ct|q| 



These relations are extensively discussed in the work of Kadanoff and Falko. pi 



IV. SOUND PROPAGATION IN THE COLLISIONLESS LIMIT 



As discussed in the introduction, in a metal like Sr2Ru04 the coUisionless limit is somewhat difficult to reach 
because we need the sound wave with the frequency u ^ 0(1) GHz. However, we will also see that this is the regime 
where one has the best chance to observe the effect of the collective mode. 

The sound velocity shift and the attenuation coefficients can be calculated by looking at the autocorrelation func- 
tions, {[Tij,Tij]), of stress tensor, r^, as discussed in the previous section. Wc will use the finite temperature Green's 
function technique ||l7|| to compute these correlation functions. The single particle Green's function, G{iujn, k), in the 
Nambu space is given by 

G^^(«a;„,k) = iujn - ■^kPs - A(fc • p)ai, (14) 

where pi and ai are Pauli matrices acting on the particle-hole and spin space respectively, ujn — {2n + l)7rT is the 
fermionic Matsubara frequency, and = k^/2m — p,. Then, for example, the irreducible correlation function can be 
computed from 

(N:Tij])oo(j'^i/,q) = 7'^^Tr[(^^^^ P3G'(p,w„)p3G'(p-q,iw„-zw^)] , (15) 

n p 

where = 2vt:T is the bosonic Matsubara frequency. 
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A. Longitudinal sound wave 



Let us consider the longitudinal wave with u || q || x, which corresponds to the Cn mode. Since the stress tensor 
couples to the density, the autocorrelation function {[hi, hi]) is renormalized by the Coulomb interaction. In the long 
wavelength limit and for s = w/wi^jql <^ 1, the renormalized correlation function {[hi, /i;])o can be reduced to 



Rc{[hi,hi])o^ (- 



Re^(cos(2</)),cos(20)) = ^A^(O) 



^-2.^(1-/) 



where 



{A,B) ^tY^Y. Tr[^P3G(p, uo,,)Bp^G{v - q, iuJn - ioJu)\ 



(16) 



(17) 



with A and B being some functions of (f> or operators. Here is the angle between p and q, N{Q) ~ m,j2-K is the 
density of states at the Fermi level and / is the superfluid density in the static limit (lo <C Wi?|q|) given by 



/ = 27rTA2 Y 



(18) 



n=0 



The derivation of the result in Eq.(jl^) is given in Appendix B 2. Therefore, the sound velocity shift, 8Ci, is given by 



(19) 



where A/ = / {STrmrriio^^Cf) is the longitudinal coupling constant. Here we set s = uj/vp[q[ — Ci/vp, where C; is 
the longitudinal sound velocity. In Sr2Ru04, s — 10^^ ^ \ which is very different from s ^ 1 of ^He. 

Now let us consider the correction due to the collective modes. The additional renormalization of {[hi, hi]) (in the 
Cii mode) due to the collective mode is computed in Appendix B 3 and the result is given by 



{[hi,hi]) = {^ 



Pf 
Am? 



i-2.2 
2 ' 



1-/- 



f{vF\ci[? 



4{(wF|q|)2+4A(r)2-2u;2} 



+ 1 —ai{uj) 



(20) 



As one can see from the above equation, there is no resonance because lo <^ wfIqI- However, we will be able to see a 
shadow of the collective mode in the sound velocity change, which we discuss in the following. 



In the limit s <C 1 and setting s = Cf /vp, the above equation leads to the sound velocity shift given by 



sell 
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/ 



A{l + i2A/vFq)^} 



(21) 



Note that the sound wave gets soften more by the collective mode. In Fig. 1, we show / = 1 — / — 4[i+(2A(T)/t;F|q|)^] 
for t;i;'|q|/A(0) = 0, 1,2,3 for 0.7 < t < 1.0 where t = T/Tc- Note that the coupling to the collective mode can be 
observed for w_F|q| ~ 2 — 3A(0) which corresponds to uj ^ 0(1) GHz. 
The attenuation coefhcient, ai, is given by 



duj'- 



(w) ~ u "~ A2^(cj'+t^)2- A^ 



f tanh 



UJ + uj' 
2T 



tanh 



-) 

2T J 



(w-2A)- 



duj'- 



L0'{uj' — Lo) 



Vt^'2- A2^(tj'-tj)2 



A2 



f tanh ■ 



2T 



(22) 



where a" is the attenuation coefficient in the normal state. This form is the same as the one in the s-wave supercon- 
ductors. 

We can carry out a parallel analysis for the longitudinal wave with u || q || z which corresponds to the C33 mode. 
Unfortunately, this sound wave does not couple to the clapping mode. Therefore, the velocity shift of this sound wave 
is simply given by 



^Of^ — \33 



1 



--2 (1-/) 



(23) 
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B. Transverse sound wave 



Here we consider first the Cge mode that has u || y and q || x. In this case, the sound velocity change can be 
obtained from the evaluation of ([/it, ht]). Assuming that the current contribution is negligible at low frequencies and 
following the same procedure used in the case of the longitudinal sound wave, we obtain 



v66 



'66 



2 \VF 



— 1-/ 



/ 



A{l + i2A/vFqn 



(24) 



where the At = Pp/ {8t. 



iCf) is the transverse coupling constant. Note that the transverse sound velocity increases 



upon entering the superconducting state. However, due to the fact that the transverse velocity is rather small and the 
coupling to the electron system is also weak compared to the longitudinal case, it will be hard to observe the change 
of the transverse sound velocity at T = Tc. 

Another transverse sound mode, C44, that has u || z and q || x, does not couple to the clapping mode. Thus the 
sound velocity change in this case is given by 



SCf^ 



v44 



1 . f c^^ 



(25) 



V. THE DIFFUSIVE LIMIT 

In the frequency range kHz ~ MHz, the diffusive limit is more realistic. In this limit, the incorporation of the 
quasi-particle damping is very important. Here we assume for simplicity that the quasi-particle scattering is due to 
impurities. Unlike the case of s-wave superconductors, we treat the impurity scattering in the unitary limit. Then 
the effect of the impurity is incorporated by changing w„ to d)„ (renormalized Matsubara frequnecy) in Eq. (^) [p^ . 
The impurity renormalized complex frequency, a)„ is determined from 



(26) 



where F is the quasi-particle scattering rate and the quasi-particle mean free path is given by / — vp/i'iT). 

In order to compare the results in the normal state and those in the superconducting state, let us first work out 
the correlation functions in the normal state, where A = 0. 



A. Normal state 



We can use Eq.(p35|) to compute {[hi, hi])o- In the limit of uj, vpQ ^ 2r, we get 
(cos(2(/)),cos(2<?!))) = (cos2(20) ( 1 - 



2ir - c 



^2r' 2r 



(27) 



One can show that (cos(2(/)), 1) is of higher order in ijj/2r and WF|q|/2r while (1, 1) « 2(cos(2(/)), cos(2(/))) to the lowest 
order. Thus, as in the previous section, {[hi, hi]) is well approximated by (cos(2(/)), cos(20)) times a multiplicative 
factor. This gives us 



1 



4r^ 



a, 



(28) 



where uj is set to C;|q|. 

It is not difficult to see that the results of the transverse sound wave is essentially the same as the longitudinal case 
up to the lowest order with a simple replacement of A/ and C; by At and Ct- Therefore, in the diffusive limit, the 
longitudinal and transverse sound velocities have the same form with different coupling constants. 
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B. Superconducting state near Tc 



Now we turn to the case of the superconducting state near Tc, where the correlation functions can be computed 
from Eg. (p6|) after replacing ijj„ by (D„. In this section, we will assume 47rTc ^ 2T vf\(\\ and use 2^ 1 near T^- 
As in the previous sections, the leading contribution in {[hi,h{\) can be computed from (cos(2(/)), cos(2(/))). 

After some algebra, we finally obtain 



(cos(20),cos(20)) 



.1. s " 
"2^2?) 



1 



1 - 



A2 



1 



27rrr 



^ ^ 2 27rT 
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2^ 



SttT ^ ^ 2 
AttT ^ ^ 2 



27rT 

r 



27rr 



(29) 



where 



7 \ n 00 ^ 



fc=0 



(z + fc)"+i 



(30) 



Here ■(/'("^(z) is the poly-Gamma function and V'(^) is the di-Gamma function. This leads to 



V2r 



= 1 



^2 27rr- 



2ttTT^ ^2 2ttT' 



(31) 



where «„ is the ultrasonic attenuation coefficient in the normal state. Here we combined the subscripts I and 
because the above analysis applies to the case of the transverse sound wave as well. Only the coupling constants Xi^t 
are different. In particular, when T /2'nT w T /2'nTc <C 1, the above equations can be further reduced to 



SCi.t . 1 



Ci^t 



1 - 



47r3 T / T 

7C(3)r V^T'c 



= 1 
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(32) 



Note that the sound velocity change and the attenuation coefficients decrease linearly in (1 — T/Tc) as T —> Tc. 
This result, Eq. (^2|), for the longitudinal sound wave is consistent with the experimental observation reported in [T^ ] 
However, in the experiment, the transverse sound velocity does not show any change across Tc- The absolute value 
of the transverse sound velocity is much smaller than the longitudinal one and the coupling to the electron system is 
also much weaker than that of the longitudinal sound waves. Therefore, it is difficult to observe any change at T — Tc 
for the transverse wave, which may explain the experimental results. 

Here we neglect the coupling to the collective mode. Indeed, even in the diffusive limit, Cn and Cgg modes do 
couple to the collective mode. However, our investigation showed that the coupling to the collective mode in these 
cases is almost impossible to detect although we do not present the details of the analysis here. 



C. Ultrasonic attenuation for all temperature regimes 



The general expression of the sound attenuation coefficient for T < Tc can be obtained by following the procedure of 
Kadanoff & Falko, and Tsuneto |16|. To obtain the ultrasonic attenuation coefficient, at, we compute the imaginary 
part of the correlation function {[Txy, T^y]). We finally arrive at 



duj — 



.g(w) 1 + yV2 - Vl + y' 



(33) 



where y : 
given by 



and nF(w) — Yjie^l^ + 1) is the Fermi distribution function. The coherence factor 5(a)) is 
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+ ' (34) 



where i = d)/ A is determined from 



u r - 1 , ^ 



Similar analysis can be also done for ai. 

In the limit of |q|Z ^ 1, the above result leads to the following ratio between the attenuation coefficients in the 
superconducting state, aj^t, and the normal state, a"^. 



Notice that the Eq. (36) applies for both of the transverse and longitudinal sound waves. This result is evaluated 
numerically and shown in Fig. 2 for several T/Tc where Fc = A(0)/2 is the critical scattering rate which drives Tc to 
zero. 



VI. CONCLUSION 



We have identified a unique collective mode called the clapping mode in a p-wave superconductor with the order 
parameter given by Eq.(l). This collective mode couples to the sound wave and affects its dynamics. 

The effect of the clapping mode on the sound waves was calculated in the coUisionless limit. However, unlike the 
case of '^He, the detection of the collective mode appears to be rather difficult. One needs, at least, the high frequency 
experiment with oj ^ 0(1) GHz. 

In the diffusive limit, we worked out the sound velocity change near T = Tc and found that it decreases linearly 
in I — T/Tc which is consistent with the experiment reported by Matsui et al ||l5|. We also obtained the ultrasonic 
attenuation coefhcient for the whole temperature range, which can be tested experimentally. On the other hand, the 
coupling of the collective mode is almost invisible in the diffusive limit. 
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APPENDIX A: 

1. Clapping mode and its coupling to the stress tensor 

The fluctuation of the order parameter corresponding to the clapping mode can be written as i5Ap3 ~ e^'^'^'^uip^,. 
The relevant correlation functions for the couplings are 

(^A, cos(2(/i))(iLjy,q) = ^ Tr[(5Ap3G(p, w„)cos(2(/))p3G'(p - q,iw„ - iuj^)] , 

n p 

((5A,<5A)(zw„q) j]Tr[,5Ap3G(p,w„)(5Ap3G(p-q,zu;„-iw,)] . (Al) 

n p 

After summing over p, we get 
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{6A,cos{2(t))) = /ttTN{0)Y^ 



2 

2 , 



where (• • •) on the right hand side of the equations represents the angle average. Now summing over w„ and analytic 
continuation iuji, lo + i6 lead to 

{SA,cos{2^))^N{Q){^F), 

(■2 , r, A 2 _ 2 

(JA,<5A) =g-i-7V(0)(^— F) , (A3) 

where F is given by 



FiujX) 4A (C W j^^2_^2)2+4£;2(^2_^2)_^4(^2A2]2_16^2£;2(^2_^2)2 ■ (^4) 

In the limit of a; ^ ^'fIqI, the contribution (in {[hi, hi]); see Appendix B 3) due to the coupling with the clapping 
mode becomes 

(<5A,cos(2<^))2 co'{Ff 



5-1 - ((5A, ,5A) ' ' 4((C2 - 2A2 - uj^)F) 

where / = limq^olhiiw^o(-P') is the superfluid density and given by Eg. (p^) . We can see that the frequency of the 
clapping mode is given by \/2A from {5 A, SA). 

APPENDIX B: 

1. Longitudinal sound wave in the collisionless limit 

The longitual sound velocity shift is given by the real part of {[hi, hi]). The irreducible correlation function for the 
stress tensor can be obtained from 

4 

{[txx,Txx])oq =T'^'^Tr[^{cos (j))'^P3G{p,ujn)p3G{p ~~ q,iu;n- iujjy)] , (Bl) 

n p 

where (p is the angle between p and q. Since the stress tensor couples to the density, the correlation function is 
renormalized as 

/\h h ]\ - l\h h n ^ V{(i){[hi,n]){[n,hi]) 

{[hi,hi\)ii - {[hi,hi\)QQ ^ , — , (B2j 

1 - V{q){[n,7i]) 

where V{q) = 27re2/|q| is the Coulomb interation. This equation can be simplified in the long wave length limit 
(|q| -> 0) as 



{[hl,hl\)QK. \-\ [{[Txx:Txx])m 77 77 • (B3j 

\Uli/ {[n,n]) 
It is useful to define the following quantity for notational convenience. 

{A,B)=tY,Y. Tr[>3G(p, c^„)Bp3G(p - q, iu,, - iuo,)] , (B4) 
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where A and B can be some functions of </> or operators. Using this notation, Eq.(B3) can be rewritten as 



([^,./>,i).»(i)'ft 

= (5) 



(cos^fli, cos^c: 



(cos^(/!>, 1) (1, cos^ 



Am? 



(cos(20),cos(20)) 



(1,1) 

(cos(2</>),l)(l,cos(2</>)) 

a;i^ 



(B5) 



Then, each correlation function can be computed from 

n p 

(1, cos(20))(ic^„ q) = r ^ Tr[cos(2(/))p3G(p, ujn)p3G{p - q, lu 

n p 

(cos(20), cos{2(j))){iuj^,q) = T ^ ^ Tr[cos^(20)p3(TiG(p, u;„)p3G(p - q, iw„ - ioji,)] 



, tUJn — lOJv 



(B6) 



n p 



Summing over p leads to 

(l,l)-(7rriV(0)^|l- 



(B7) 



where C, = v^? • q and A^(0) — mjlix is the two dimensional density of states. Similar equations can obtained for 
(cos(20), 1) and (cos(2(/)), cos(2(/))) with additional angle factors, cos(2(/)) and cos2(20) respectively. After summing 
over and analytic continuation iw^ — > + ?(5, we get the following results in the limit of w ^ I'Flql- 



(1,1)=7V(0)( 



(cos(2(?!)),l) = JV(O)(cos(20) 



iV(0) 



(cos(20),cos(20)) = 7V(0)(cos2(2^) 



iV(0) 
«iV(0) 



2\/r~~ 



(B8) 



We find that the second term in the last line of Eq.(B5) is of higher order in u}/vp\c^{= s) so that we can ignore 
it. In Sr2Ru04 or metals, s ^ 1. Thus the effect of the coupling to the density is merely to change the vertex 

2 2 

associated with Txx from ^cos^4> to |^cos(20) as far as the lowest order contribution is concerned. Evaluation of 
(cos(20), cos(20)) leads to 



Ke{[hi,hi\)a^ ( - ) Re 



Am? 



(cos(20),cos(20)) = - 



Pf 



Am? 



7V(0) 



^-2.^(1-/) 



(B9) 



where / is the superfiuid density. 



2. Contribution coming from the coupling to the clapping mode 



The correction due to the collective mode leads to the renormalized correlation function as follows. 



{[huh]) = {[huhi])o + 



- ([<5Ap3,<5Ap3 



(BIO) 



where g is the coupling constant between the stress tensor and the collective mode. (5Ap3 represents the fluctuation 
associated with the clapping mode. Using the fact that (l,e2*"^) — and (5A ~ e^^'^ax, the above equation can be 
further reduced to 



huh,]) = 



2 P% 



2 

Pf 
4m2 



cos(20),cos(2(/))) 

2 

vf 



(cos(20),(5A)(5A,cos(20)) 



2 



1-/ 



(,5A, (5A) 



4{(z;f|q|)2+4A(T)2-2a;2} 



+ ,!!!i2ii^aiM , (Bii) 
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where / is the superfluid density and given by Eq.(pJ 
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FIG. 2. The normalized attenuation coefficient as a function of the 
reduced temperature t = T/T^ for T/Fc = 0, 0.1, 0.2, 0.3, 0.4. 
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